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Double-dipole excitations in 208 Pb are analyzed within a microscopic model explicitly treating 2p2h- 
excitations. Collective states built from such 2p2h-excitations are shown to appear at about twice 
the energy of the isovector giant dipole resonance, in agreement with the experimental findings. The 
calculated cross section for Coulomb excitation at relativistic energies cannot explain simultaneously 
the measured single-dipole and double-dipole cross sections, however. 



1 Introduction 

Double giant dipole resonances (DGDR) have been observed in pion double-charge- 
exchange reactions at LAMPF as well as in high-energy heavy-ion collisions at GSI 
P] . In the latter case there is a strong focusing of the electromagnetic field in the target 
rest frame. This greatly enhances the field intensity in the vicinity of the target nucleus, 
thus increasing the probability for two-photon absorption from the ground state. 

The global parameters of the DGDR, i.e. the resonance energy Ejjqdr, the width 
Fdgdr and the cross section a can be summarized as follows ||: 1) The resonance 
energy is about twice as large as that of the isovector giant dipole resonance (GDR), 
Edgdr — 2Egdr- which suggests that can be interpreted as an independent two-phonon 
state. 2) The observed values of the width Tdgdr are bracketed by V2Tqdr and 2Tgdr- 
Also this feature supports a harmonic picture. 3) The measured cross section <7 exp is 
larger than the theoretical estimate <7th- The values of the ratio <7 eX p/<7th are scattered 
between 1 and 5 depending on the nuclei considered and on the theoretical analyses. 
This puts into some doubt to interpret the DGDR as an independent two-phonon state. 

To reach a quantitative understanding of the DGDR characteristics we shall carry out 
a microscopic calculation of the DGDR in 208 Pb. The structure of double-dipole states 
are investigated by the diagonalization of a model Hamiltonian in the space of lplh- 
and 2p2h-excitations. For comparison with experiment the cross section for Coulomb 
excitation, based on the second-order perturbation theory, will be estimated. 



2 Double-Dipole Strength Distribution in Pb 
2.1 Microscopic Model 

The nuclear structure model has been developed in our previous work j| on double-dipole 
excitations in 40 Ca. We give a brief outline stressing some new aspects in the application 
to heavy nuclei such as 208 Pb. The eigenstates |n) of the nuclear Hamiltonian H$ are 
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determined by diagonalizing Hq in the space of lplh- and 2p2h-excitations. As a model 
space of single-particle states, we take into account three major shells of the Woods- 
Saxon potential on both sides of the Fermi level. The continuum states are treated by 
discretization through an expansion in a harmonic oscillator basis. The number of 2p2h 
states constructed with these single-particle states exceeds several hundreds of thousand. 
Most of them, however, have vanishingly small matrix elements for the double-dipole 
transition from the ground state. Thus we select only those 2p2h-states whose energy 
is below 40 MeV, and which exceed a lower limit in the magnitude of the double-dipole 
matrix element, 

\{2p2h\Db\Q)\ 2 /Y, l(2|^D|0>| 2 > 5 x 1(T 5 . (1) 



Here, the isovector dipole operator, b, is defined as 

- N* Z" 

i=l i=l 

With these conditions, we can reduce the number of 2p2h-states to 996 for the J 71 " = + 
double-dipole transition and to 2011 for the 2 + transition. 

The mean energy and total strength of the GDR in 208 Pb, calculated in the usual 
lplh-RPA with the nuclear Hamiltonian jl, including a density-dependent zero-range 
interaction, are 10.3 MeV and 98% of the TRK sum-rule value, respectively. Both are 
smaller than photo-neutron data || which yield 13.4 MeV and 134%. Therefore we 
enhance the theoretical results by a scaling of the single-particle energies according to 
£sp = £ws/( m * / m ) with the effective mass m*/m = 0.75. By this procedure, we obtain 
the mean energy and the total strength as 13.5MeV and 131% of the TRK value, in good 
agreement with experiment. 



2.2 Strength Distribution 

The strength distributions of the double-dipole transition are shown in Fig. |l]. Figs. 1(a) 
and 1(b) (1(c) and 1(d)) represent the unperturbed and the perturbed strength dis- 
tributions for the + (2 + ) double-dipole transitions, respectively. Both for + and 2 + 
double-dipole transitions, the perturbed strength concentrates in the region of 25-30 
MeV. There appear several states which carry about ten times more strength than the 
strongest unperturbed ones. 

The mean energy, (E), and the width, T, are estimated in terms of the energy 
moments 




dEE k S DD (E), S DD = J2\(n\bb\0)\ 2 5(E-E n ), (3) 



as 

{E) DGDR = m 1 DD /m {i DD and T DGDR = ^Jm 2 DD m° DD - (m l DD ) 2 /m° DD , (4) 

and similarly for the isovector dipole excitations. The (E) dgdr and ^dgdr are 26.96 
MeV (26.53 MeV) and 3.48 MeV (3.36 MeV) for 0+ DGDR (2+ DGDR), respectively. 
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Figure 1: Double-Dipolc Strength Distributions, (a) Unperturbed + DD strength, (b) Perturbed 
0+ DD strength, (c) Unperturbed 2+ DD strength, (d) Perturbed 2+ DD strength. 

The mean energy of the DGDR is about twice as large as that of the GDR. The calculated 
width, which corresponds to the fragmentation width (Landau damping) is larger than 
that of the GDR (T GDR = 2.70 MeV) by the factor of 1.24-1.29. We note that the 
observed width includes, in addition, the escaping and spreading width which are not 
included in the present calculation. 

3 Cross Section for the Coulomb Excitation 

3.1 Second- order Coulomb -Excitation Amplitude 

In the semiclassical approach [^,@], the cross section for Coulomb excitation is given by 

fOO 1 

0W/ = 2tt / bdb—— J2 M 2 - ( 5 ) 

The projectile is assumed to move on a straightdine trajectory with impact parameter 
b. The minimum impact parameter R m in is a key parameter, since the cross section 
sensitively depends on it as will be shown later. 
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The main components of the DGDR are 2p2h-states, which are excited via a two- 
body operator of double-dipole character (DD) from the ground state. For Coulomb 
excitation the precise form of the excitation operator is derived from the 2nd-order 
amplitude 

a (2) = _L_ (/ | / dt / dt > ^iHot/Ky^-iH.t/n^iH.t'/Hy^y-iH.t'/n^^ (g) 

(in) J — oo J— oo 

where Hq and V are the nuclear Hamiltonian and the Coulomb interaction between 
the projectile and target nuclei, respectively. To obtain the two-body operator, we 
exchange the order of V(t) (V(t')) and exp (+iHQt' /h) (exp (— iHot/h)). This yields a 
series expansion in the times t and t' which includes single- and multiple-commutators 
of Hq and V . In the limit of a 'fast collision' Q, we need only to take into account the 
lowest order in time. Finally the expression for the second-order amplitude of the double 
EX' transition is given by 

(2) Ze. 2 ^ (-)" 

+ (/|X:^'(A,-M)K)f V f^^U^'^-q^ + q)}, (7) 

with 



T ^W 2( ^ 2) = [S„ lXl (u; 1 )®S W2X2 (u; 2 )} Xlji (8) 
E^i Wa (wi,w 2 ) = PW"i)®** a A>2)]^ (9) 

tf£%(\u) = [MfrMsMi^)]^ (10) 

/C^ 2 2 (A, M ) = [[HcMCTnAOlsMCTraAa) 

-M(7riAi)®[F ,^(7r 2 A 2 )]] AM . (11) 

Here, .M(7rA) denotes electric (ir = E) or magnetic (tt = M) multipole operators while 
^ a and R^x are orbital integrals || . The principal value integral in Eq. [7| results from 
a step function, which appears when the upper limit of the integral in Eq. || is extended 
to infinity. In the present estimate of the cross section, we neglect this principal value 
integral. 

3.2 Cross Section 

The calculated cross section for Coulomb excitation of the DGDR with a 208 Pb projectile 
incident on a Pb target at 640 MeV/A is plotted as the solid line in Fig. ||. The dashed 
and dotted lines display the cross sections for the GDR and the GQR (giant quadrupole 
resonance), respectively. For comparison the former has been multiplied by the factor 
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of 0.1. To account for experimental energy resolution, the calculated cross sections have 
been smoothed by using a Breit-Wigner function of width 1.5 MeV. The main peak of 
the DGDR appears in the region of 25 - 30 MeV, just above the broad isovector GQR 
around 20 MeV. The peak energy of the DGDR is about twice that of the GDR. 
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Figure 2: Cross Section for the Coulomb Excitation. The solid line denotes the cross section of 
the DGDR while the dashed line displays the GDR cross section (multiplied by 0.1). In addition, 
the cross section for the GQR is shown as the dotted line. The minimum impact parameter is 
the value of case (b) (see text). 



Next we discuss the dependence of the energy-integrated cross section on the mini- 
mum impact-parameter, R m i n . We parameterize R m in as R m in = r o{A l J 3 + A l J 3 ), where 
At and A p denote the mass number of target and projectile nuclei, respectively. Three 
cases are considered : (a) ro = 1.50 fm, (b) r$ = 1.31 fm and (c) ro = 1.20 fm. The 
values for (a) and (c) have been used previously by Ponomarev et al. |lG| in the study 
of the DGDR in 136 Xe, while case (b) represents the choice of Boretzky et al. || (Fig. Q 
corresponds to case (b)). The cross sections of the GDR and the DGDR in 208 Pb inci- 
dent on various nuclei are given in Table |l] where the three cases for R m in are compared 
with the measured cross sections ||. The intermediate value of R m i n (case (b)) repro- 
duces the measured cross sections of the DGDR fairly well, but overestimates those of 
the GDR. On the other hand, the larger value (case (a)) reproduces the measured cross 
sections of GDR, while underestimating those of the DGDR by the factor 2-3. In our 
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calculation it seems impossible to simultaneously explain the measured GDR and DGDR 
cross sections with the same value of R m in ■ 

In the particle-vibration coupling calculation of the DGDR in 136 Xe by Ponomarev et 
al. [1C], similar results for the cross section were obtained, but the discrepancy between 
the measured cross section and their estimate is larger than that in our results. On the 
other hand, the simple folding model analysis of the DGDR in 208 Pb by Boretzky et al. 
||] predicts only a 33% enhancement of the DGDR cross section while reproducing the 
GDR cross section. It remains an open question whether this difference between results 
from the microscopic models and the folding model originates from the nuclear structure 
or from the treatment of the reaction mechanism. 



Table 1: Cross section for Coulomb excitation of the GDR and DGDR in a 208 Pb projectile 
incident on U, Pb, Ho and Sn targets at 640 MeV/A. 



Cross 
Section (b) 


U 


Pb 


Ho 


Sn 


GDR 


DGDR 


GDR 


DGDR 


GDR 


DGRD 


GDR 


DGDR 


(a) 


4.18 


0.222 


3.47 


0.162 


2.49 


0.091 


1.52 


0.038 


(b) 


5.34 


0.500 


4.42 


0.361 


3.16 


0.201 


1.91 


0.083 


(c) 


6.22 


0.827 


5.12 


0.590 


3.63 


0.321 


2.18 


0.129 


exp 


3.66 


0.51 


3.28 


0.38 


2.47 


0.28 


1.45 


0.07 



4 Summary 

We have carried out a microscopic calculation of double-dipole excitations in 208 Pb in a 
realistic model space of lplh- and 2p2h-states. Since the number of 2p2h-states in heavy 
nuclei such as 208 Pb, is prohibitively large we have introduced the selection procedure 
given by Eq. [l] as well as a truncation in the 2p2h energy. Although the width of DGDR 
may not be totally reliable because of this procedure, we believe that both the mean 
energy and the integrated strength are realistic. 

The double-dipole strength is shown to be concentrated in an energy region twice 
that of the isovector giant dipole resonance thus indicating that anharmonicity effects 
are quite small, in good agreement with experiment. 

The Coulomb excitation cross section has been calculated based on the second-order 
perturbation. It depends sensitively on the minimum impact parameter, which appears 
in the semiclassical treatment of the cross section for Coulomb excitation. We have 
estimated the cross section of the GDR and DGDR with three choices for the minimum 
impact parameter, available from the literature. None of these choices can explain the 
measured cross sections simultaneously and hence the discrepancy between the measured 
cross section and theoretical estimates, previously observed, remains. 
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